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X 1 Abstract 

«?; 

Uncertainties in the NLO calculation of the inclusive jet cross section due 
to the choice of renormalization scale, parton distribution functions and clus- 
tering algorithm are explored. These are found to be similar in size to the 
current experimental uncertainties of the measured inclusive jet cross section 
at D0 and CDF. 
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I. INTRODUCTION 



Measurement of the inclusive jet cross section in proton anti-proton (pp) collisions con- 
stitutes a strong test of the pdictions of perturbative quantum chromodynamics (QCD). 
Deviations of the theoretical cross section from the experimentally observed cross section 
may be evidence for physics beyond the Standard Model. In particular, the presence of 
quark compositeness would enhance the cross section at high values of transverse energy 
(E T ). 

Recently, the inclusive jet cross section has been measured by the CDF JTJ and D0 |2| 
Collaborations with systematic uncertainties ranging from 10% to 40% as function of E^. 
With the improvement in the experimental accuracy of this measurement, it is worth inves- 
tigating the accuracy of the next to leading order (NLO) QCD predictions 0-01 . 

In this paper the uncertainties in the NLO QCD inclusive jet calculations are explored 
using two available programs: Jetrad Q a complete 0(a s 3 ) event generator, and EKS || 
a complete 0(a s 3 ) analytical calculation of the inclusive jet cross section. Both programs 
require the selection of a renormalization and factorization scale (typically chosen to have the 
same value, /i), a set of parton distribution functions (PDF), and a jet clustering algorithm. 

II. DISCUSSION OF THE THEORETICAL PREDICTIONS 

If QCD could be calculated to all orders the results would be independent of the choice 
of renormalization scale. Because the inclusive jet cross section has been calculated only 
to NLO (0(a s 3 )), the choice of renormalization scale does affect the result. The authors 
of Jetrad have provided several choices for the renormalization scale, we have chosen to 
investigate a scale proportional to the E T of the leading jet after parton clustering (// = 
AE™ X , where A is a constant typically chosen to lie in the range 0.25 < A < 2). The 
authors of EKS prefer an alternative definition of the renormalization scale: the Ex of 
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each jet in the event (/i = A-E^Q- An alternative scheme (available in Jetrad) uses the 
center-of-mass energy of the two outgoing partons as the renormalization scale (/i = C\fl = 
C^xix 2 s where C is a constant chosen to lie in the range 0.25 < C < 1, Xx — Yl E Ti e Vi /y/s, 
X2 = J2 ETie~ Vi I \/s, i = 1 . . . n where n is the number of jets in the event, rj = — ln[tan(#/2)] 
and 9 is the polar angle relative to the proton beam )Q. The effect of these scale choices on 
the inclusive jet cross section is discussed. 

The standard Snowmass clustering algorithm Q combines two partons into a single 



jet if they are both within 1Z = y r/ 2 + <p 2 = 0.7 of their weighted center (where is 
the azimuthal angle). An additional constraint on the parton clustering is applied which 
requires that the two partons be closer than lZ scp = 1.3x71 f|f]. The value of 7Z sep has been 
chosen to match the experimental jet splitting/merging parameters used in the jet clustering 
algorithms ||. The effect of using 7Z sep will be discussed. 

The uncertainty in the calculation of the inclusive jet cross section resulting from parton 
distribution functions will be divided into three parts. 

The first is due to the choice of the PDF family. This choice is associated with the 
selection of data used to determine the PDF and the functional form used in the fits. The 
variation of the inclusive jet cross section has been studied using a selection of modern PDFs: 
cteq3m ||, cteq4m, cteq4hj |10[ and mrsa' fill . 

The second category of PDF uncertainties results from the value of the strong coupling 
constant (a s ) used in the PDF. Usually a free parameter in the PDF fit, a s can be fixed 
to a pre-determined value. The effect of a s variation on the PDF is examined by using 



l K version of EKS that uses the renormalization scale fi = AE™ ax is also available. 

2 The choice of renormalization scale [i = C\pk is somewhat unnatural for the inclusive jet cross 
section which is dominated by i-channel exchange. It has been included to study the effect of an 
extreme choice of scale and for comparison with previous two-jet mass analyses ||. 

3 The Snowmass algorithm corresponds to TZ sep = 2.0 x 1Z 
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the cteq4a series in which a s is fixed to values ranging from 0.110 to 0.122 at Mz- In 
comparison, a s (M z ) = 0.116 for cteq4m. 

Finally, Jetrad and EKS use different strategies to evolve the PDF in x and Q 2 (where 
x is the momentum fraction carried by the parton and Q is the characteristic energy scale of 
the process, typically chosen to be the momentum transfer). The Jetrad program uses the 
strategy as implemented by the MRS group for evolving all PDFs (Jetrad also uses 
CTEQ PDFs generated using the CTEQ evolution package). EKS implements the PDFs 
by interpolating from a table of values that were generated directly from the original PDFs. 
While these implementations are theoretically identical, small differences can be produced 
by the numerical accuracy of the program. 

The uncertainties will be determined at yfs = 1.8 TeV and 630 GeV and compared to 
a reference model. For this study, the reference model will be the Jetrad calculation for 
the pseudorapidity range \rj\ <0.5, /i = 0.5E™ ax , 1Z scp = 1.3 and the CTEQ3M PDF evolved 
using the CTEQ method (see Fig. [I]). The comparisons with other theories will be given by: 

/ Theory \ 
y Reference Theory J 

This ratio is fitted with a third degree polynomial yielding a smooth curve. In most cases 
the ratio can be fitted with a resulting \ 2 P er degree of freedom less than one. 



III. INCLUSIVE JET CROSS SECTION AT y/S = 1.8 TEV 

The difference between using the standard Snowmass clustering algorithm and the mod- 
ified algorithm with TZ sep = 1.3 is shown in Fig. |2|. The effect ranges from 8% at 50 GeV 
decreasing to 5% at 500 GeV. Because the value of 7Z sep is selected to reflect the experimen- 
tal clustering algorithm, the uncertainty resulting from its use is much smaller than 5%. A 
more appropriate variation of the value of TZ sep is from 1.2-1.4 resulting in variations in the 
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cross section of less than 1% 

The effect on the cross section due to the choice of the renormalization scale was studied 
by using Jetrad with several different values of fi. First, the cross section was calculated 
using n = 0.25, 0.75, 1.0 and 2.0 E™ x . These are compared to the cross section with fi = 0.5 
^max j n pjg^ g rj^ e crosg sec ^i on i s largest^] for fi = 0.5E™ ax and is reduced by 5-10% with 
some Et dependence for fi = 1.0E™ ax and 0.75E™ ax . The cross section is approximately 
10% below the reference model for /i = 0.25i?™ ax and 15-20% below the reference model for 
fi = 2.0E^ ax with some dependence. 

Figure |] compares the predictions for alternative choices of fi, and V§. The choice 
/i = O.SE^ (calculated with EKS0) is compared to fi = 0.5_E™ ax and shows a 5% difference 
at an E^ of 50 GeV dropping to less than 1% at 500 GeV. Also shown in Fig. |] is a 
comparison of the Jetrad calculations for fi = 0.25, 0.5, 1.0 compared with /i = 0.5 
Ef ax . The effect is approximately 20-25% at 50 GeV decreasing to 10% at 500 GeV with a 
strong dependence on the choice of scale used. In summary, there are only small differences 
between the scale choices of E^ and i?™ ax but large differences when the scale is changed 
from E"c to V§. 

The effect of different PDF choices is depicted in Fig. [5[ The cross section calculated 
with cteq4m is a few percent lower than the cross section calculated using cteq3m. These 
differences result from a change in the parameterization used to model the gluon distributionP] 

4 The uncertainty resulting from the choice of !Z S ep will not be considered in the remainder of this 
note. 

5 The inclusive jet cross section at fi = 0.5 £™ ax is a maximum since this is a point of minimum 
sensitivity for the calculation, see |l^] for a discussion. 

6 Jetrad does not implement fi = AE^ . 

7 CTEQ3m uses a more restrictive gluon parameterization: G(x,Qo) = Aqx Ai (1 — x) A2 (l + A%x) 
while cteq4m uses: G{x,Q ) = A x Al (l - x) A * (1 + A 3 x Ai ) 
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of the cteq4m PDF and the inclusion of additional data sets. These data sets include more 
precise deep inelastic scattering (DIS) data from NMC [13j and HERA ||I~^ , |T5|| , the inclusive 
jet cross section measured by CDF fiyiBfl and the preliminary inclusive jet cross section 
measured by D0 |I7j]. These changes lead to a change in the optimal value of a s from 0.112 
to 0.116 The cteq4hj PDF, which emphasizes recent Tevatron jet data to constrain 
the gluon distribution, shows a decrease in the cross section at low Et of approximately 
5% and an increase in the cross section at 500 GeV of approximately 25%. The final 
comparison is made using mrsa', which uses a slightly different parameterization and input 
data to CTEQ3M. This results in cross sections that are similar at high and low Et and 
approximately 5% higher at 300 GeV. These differences are caused by the gluon distributions 
of the PDFs (see Fig. |6|). Note that the variations in the cross section calculations due to 
the choice of PDF are limited by the similarity of the parameterizations used to model the 
gluon distributions in the PDFs (which are not well constrained by experiment). The choice 
of an alternative parameterization of the gluon PDF could lead to larger uncertainties. 

The cross section change due to the variation of a s in the PDF is shown in Fig. [7]. 
Approximately ±5% changes are seen at low E^ which diminish as the E^ increases. The 
uncertainty in the PDF due to the choice of a s is significantly smaller than the uncertainty 
due to the gluon distributions. 

A comparison can be made between the EKS and Jetrad calculations for | rj | <0.5 using 
H = 0.5E™ aX ) = 1-3 and the CTEQ3M PDF. As shown in Fig. g, the two calculations 
differ at a level of 2-3% with some dependence on the E^. This variation is due to the 
different evolution in x and Q 2 used by the two programs. Figure |9| shows the variation 
in the Jetrad predictions due to the different evolution methods. These differences lead 
to 5% difference in the cross section with some E^ dependence, similar to the differences 
between the EKS and Jetrad programs. 



6 



IV. INCLUSIVE JET CROSS SECTION AT y/S = 630 GEV 



The studies described in the previous section were repeated for y/s = 630 GeV. The 



comparison between Jetrad and EKS with fi = AE™* is given in Fig. [H| EKS produces a 
cross section that is 5% lower than the Jetrad cross section at 20 GeV and 10% higher at 
150 GeV. These uncertainties cannot be fully explained by the difference in choice of PDF 
evolution. 

The variations in the NLO calculations of the cross section due to the choice of renor- 
malization scale and PDF are given in Fig. These variations are slightly larger than 
those observed at yfs = 1.8 TeV. 



V. THE RATIO OF THE INCLUSIVE JET CROSS SECTIONS AT \/S = 1.8 TEV 

AND 630 GEV 

Theoretical uncertainties in the NLO QCD predictions should be reduced in the ratios 
of the inclusive jet cross sections at \/s = 630 GeV and 1.8 TeV as a function of the 
dimensionless quantity X T = 2E T / y/s. 

The variations in the ratio of the cross sections due to the choice of renormalization 
scale is approximately 15% with some dependence on Xt (see Fig. [12] (a) and (b)). The 
uncertainty due to the choice of PDF is only a few percent (Fig. [T2| (c)) which compares 
to an uncertainty of up to 25% in the individual cross sections. The difference due to the 
variation of a s is not reduced by measuring the ratio of the cross sections and is still at the 
5% level (Fig. 0(d)). 



VI. CONCLUSION 

The inclusive jet cross section, predicted using the available NLO programs, has signifi- 
cant uncertainties due to the choice of renormalization scale and PDF. The overall variation 
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in the cross section can be as large as ±30%. Except for alternative implementation of PDFs 
and the evolution strategies used, Jetrad and EKS appear to be identical. 

The ratio of the inclusive jet cross sections at y/s = 630 GeV and 1.8 TeV as a function 
of Xt has uncertainties of approximately 10-20%, which is much smaller than the variation 
of the cross sections. 

Before the inclusive jet cross section can be used to test QCD or search for New Phe- 
nomenal, the theoretical predictions must improve. Most feasibly through improved mea- 
surement of the gluon distributions. 

Motivation for this work work arose as a result of discussions within the D0 QCD group 
and in discussions between that group, W. Giele, D. Soper, and the CDF QCD group. 
We thank the staffs at Fermilab and collaborating institutions for their contributions to 
this work, and acknowledge support from the Department of Energy and National Science 
Foundation (U.S.A.). 



8 This result does not not eliminate the possibility of observing new physics that produce effects 
on the inclusive jet cross section that are significantly larger than the theoretical uncertainties. 
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FIG. 1. The reference NLO calculation of the inclusive jet cross section at y/s = 1.8 TeV. 
Jetrad calculation for the pseudorapidity range \rj\ <0.5, fi = 0.5£'™ ax , lZ scp = 1.3 and the 
CTEQ3M PDF evolved using the CTEQ method. 
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FIG. 2. A comparison between the Jetrad calculations of the inclusive jet cross section with 
K sep = 2.0 and 1.3 at s/s = 1.8 TeV. 
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FIG. 3. A comparison between the Jetrad calculations of the inclusive jet cross section with 
H = 0.25, 0.75, 1.0 and 2.0£^ ax compared with n = 0.5E™ ax at yfa = 1.8 TeV. 
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FIG. 4. A comparison between the EKS calculation of the inclusive jet cross section with 
H = O.bE 3 ^ and (i = 0.5E™ ax . Also shown is a comparison between the Jetrad calculations of 
the inclusive jet cross section with fi = 0.25, 0.5 and 1.0\/I compared with \i = 0.5E™ X at yfs = 
1.8 TeV. 
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FIG. 5. A comparison between the Jetrad calculations of the inclusive jet cross section with 
the cteq4m, cteq4hj and the MRS a' PDF compared with the calculation using cteq3m at y/s 
= 1.8 TeV. 
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FIG. 6. A comparison of the PDFs used in this analysis, (a) The CTEQ3M PDF for Q = 100 
GeV (the momentum transfer) as a function of the momentum fraction carried by the parton [x). 
(b) A comparison between the gluon distributions from the cteq4m, cteq4hj and the MRS a' 
PDFs compared with the CTEQ3M PDF. (c) A comparison between the up quark distributions 
from the cteq4m, cteq4hj and the mrsa' PDFs compared with the CTEQ3M PDF. 
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FIG. 7. A comparison between the Jetrad calculations of the inclusive jet cross section with 
the cteq4a series of PDFs compared with the calculation using cteq4m at y/s = 1.8 TeV. 
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FIG. 8. A comparison between the Jetrad and EKS calculations of the inclusive jet cross 
section at y/s = 1.8 TeV. The theory parameters are \rj\ <0.5, fi = 0.5i?™ ax , lZ scp = 1.3 and the 
CTEQ3m PDF. The fit is to a third degree polynomial. 
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FIG. 9. A comparison between Jetrad calculation of the inclusive jet cross section with the 
CTEQ3M PDF using the MRS and the CTEQ evolution packages. 
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FIG. 10. A comparison between the Jetrad and EKS calculations of the inclusive jet cross 
section at yfs = 630 GeV. The theory parameters are \ rj\ <0.5, fi = 0.5i?™ ax , 1Z scp = 1.3 and the 
cteq3m PDF. 
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FIG. 11. The deviations of the inclusive jet cross section at sj~s = 630 GeV. (a) A comparison 
between the Jetrad calculations with /j, = 0.25, 0.75, 1.0 and 2.0£ , ™ ax compared with \x = 0.5^™ ax . 
(b) A comparison between the EKS calculation with \i = Q.hE^ and // = 0.5i?™ ax . Also shown 
is a comparison between the Jetrad calculations with fi = 0.25, 0.5 and l.0V§ compared with 
fj, = OSE™ 3 *. (c) A comparison between the Jetrad calculations with the cteq4m, cteq4hj 
and the mrsa' PDF compared with the calculation using cteq3m. (d) A comparison between 
the Jetrad calculations with the cteq4a series of PDFs compared with the calculation using 

CTEQ4M. 
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FIG. 12. The deviations of the ratio of inclusive jet cross sections at ^/s = 1800 GeV and 630 
GeV. (a) A comparison between the Jetrad calculations with /j, = 0.25, 0.75, 1.0 and 2.01?™ ax 
compared with \x = 0.5i?™ ax . (b) A comparison between the EKS calculation with fi = O^E 3 ^ and 
fi = 0.5£?p ax . Also shown is a comparison between the Jetrad calculations with /j, = 0.25, 0.5 
and 1.0\/j compared with = 0.5E™ ax . (c) A comparison between the Jetrad calculations with 
the cteq4m, cteq4hj and the mrsa' PDF compared with the calculation using cteq3m. (d) A 
comparison between the Jetrad calculations with the cteq4a series of PDFs compared with the 
calculation using cteq4m. 
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